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Tracing spiral density waves in M81. 
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ABSTRACT 

We use SPITZER IRAC 3.6 and 4.5/im near infrared data from the Spitzer Infrared 
Neai-by Galaxies Survey (SINGS), optical B, V and / and 2MASS Kg band data to produce 
mass surface density maps of M8 1 . The IRAC 3.6 and 4.5/im data, whilst dominated by emis- 
sion from old stellar populations, is corrected for small-scale contamination by young stars 
and PAH emission. The / band data are used to produce a mass surface density map by a 
B-V colour-correction, following the method of Bell and de Jong. We fit a bulge and expo- 
nential disc to each mass map, and subtract these components to reveal the non-axisymmetric 
mass surface density. From the residual mass maps we are able to extract the amplitude and 
phase of the density wave, using azimuthal profiles. The response of the gas is observed via 
dust emission in the 8/im IRAC band, allowing a comparison between the phase of the stellar 
density wave and gas shock. The relationship between this angular offset and radius suggests 
that the spiral structure is reasonably long lived and allows the position of corotation to be 
determined. 
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1 INTRODUCTION. 

Although many theories have been proposed, the origins and driv- 
ing forces of spiral structure in galaxies are still not particularly 
well understood. It is widely accepted that observed grand design 
spiral arms are tracers of underlying mass variations in galaxies. 
It is possible that these grand design spirals are long-lived struc- 
tures, independently able to maintain their sha pe for many galac- 
tic rotations (QSSS - iLin & Shul J1964 1 19661) ). although the ex- 
istence of isolated long-lived density waves has been questioned, 
as simula tions are only able to reprod uce transient spiral struc- 
ture (e.g. ISellwood & Carlberd ( Il984h ). If grand design spirals 
are driven rather than spontaneous, the driving force could be 
a central bar or ex ternal interactions with companion galaxies 
iKormendv & Norma n ( 1979); Elmegreen et al. ( 1982). 

It is clear from studying spiral galaxies that star formation 
can be closely linked with spiral structure, and is preferentially 
located on spiral arms. It has even been argued that the ob- 
served spiral arms are in fact an effect, rather than cause, of the 
star fo r mation (self-propagating star formation iMueller & Arnetll 
( Il976[ ): iGerola & Seided <1978r) : a chain of supernova shocks, or 
some other star-forming process, induces star formation which is 
then sheared into spiral arms by the differential rotation of the 
galaxy). This theory has been disproved for many large grand- 
design galaxies through the o bservation of unde r lying spiral arms 
in the old stellar population JSchweizea <1976h . lElmegreen et al.l 
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( ll989h . lElmegreen & Elmegreenl jl990l) . lRix & Riekd h993h ), but 
is still an intriguing possibility for flocculent galaxies, although 
some flocculent galaxies have also been shown to have weak un - 
derlying spiral structure lElmegreen et al.l <1999l) : lThornlevl ( Il996r) . 

Many surveys have been carried out using relatively large 
samples of spiral galaxies observed in th e near infrared, includin g 
work by lEl megreen & Elmegreen ( 1984), Rix & Zaritskvl ( Il995h . 
ISeigar & lames. (1998 ), and Gr osb0l et al.. C2004h which identify 
underlying density waves and place constraints on the relative am- 
plitude of the stellar arms (as a fraction of the axisymmetric com- 
ponents). The amplitude of the spiral wave is an important input pa- 
rameter in predicting the response of the gas to the stellar potential. 
In some cases, these studies also contrast the infrared stellar mass 
surface density variation with s ta r formation via optica l imag es 
JEhnegreen & Elmegreenl JI984|) : ISeigar & JamesI h998l, 12003) ). 
Other morphological studies have looked at the phases of spiral 
arms - e.g. studies of pitc h angle ( lDanveilll942l : lKennicuttlll98ll : 
iKennicutt & Hodgdll982h . Some studies find evidence of angular 
offsets between the stellar waves and star forming features or dust 
lanes has been f ound, e.g. ISeigar & JamesI ( 119981) and tentatively 
lKenmcuttlfl98lh . 

The relative locations of the spiral density maxima (and hence 
potential minima) and the shock front in the gas can be used as an 
indicator of the long-term behaviour of the stellar spiral because the 
response of the gaseous disc to the spiral potential is strongly influ- 
enced by the lifetime of the spiral pattern. There are three possibil- 
ities; i) a transient spiral where the lifetime is of order the dynam- 
ical time, such as generated spontaneously in an isolated galaxy. 
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will indu ce shocks in the gas th a t are l o cated on the poten t ial mi n- 
ima (e.g. ISellwood & Carlberd ( Il984l) : Iciarke & Gittin3 ( l2006h ). 
ii) If the spiral structure is long-lived, steady-state behaviour has 
time to emerge (a true density wave in the QSSS interpretation). 
The steady state response of the gas is to form a shock at an angu- 
lar offset which lies upstream from the density wave a t large radii 
l lRoberts|[T969l : IShu et alj[T973 : lGittins & Clarkell2004h . iii) Alter- 
natively, if - as is likely in the case of M81 - the spiral structure is 
induced by interaction with a companion, the lifetime of the pat- 
tern is likely to be intermediate between the cases i) and ii) above. 
The consequences of this scenario for the relationship between the 
arms and shocks has not been explored to date. Thus the behav- 
ior of the shock indicates whether the spiral pattern is a long-lived 
or transient feature, and if the former, potentially provides a way 
of constraining corotation without resorting to arguments based di- 
rectly on the morphology. 

M81 has of course been stu died in detail in the past, perhaps 
most notably by Ivissen jl980al lbl) through the use of HI dynam- 
ics combined wit h the density wave amplitude as determined by 
ISchweizeii 1 119761) . However, despite the extensive past research 
into spiral structure (and M81 in particular), only a small fraction 
of which is discussed above, the subject is by no means fully under- 
stood. The increased depth and resolution of infrared images pro- 
vided by the Spitzer Space Telescope offers an ideal opportunity to 
return to the subject. In this paper we produce maps of stellar mass 
surface density in order to determine the morphology of the stel- 
lar spiral structure, including the pitch angle and amplitude of the 
spiral arms. We also identify the shock front in the gas, and quan- 
tify the offset between the peak of the spiral arms and gas shock in 
order to determine the radius of corotation. 

In addition to the direct application in determining spiral struc- 
ture, a reliable and simple method of reproducing the underlying 
mass distribution in a galaxy would be very useful. Potential appli- 
cations include studies of secular evolution of galaxy discs (through 
quantifying torques), identifying weak spiral arms in optically floc- 
culent galaxies, and studying the relationship between star forma- 
tion and spiral arm strength, in addition to providing accurate maps 
of the potential distribution for input into simulations. 



2 METHODS. 

2.1 Overview of Data. 

In order to study stellar spiral structure it is necessary to produce 
maps of stellar mass surface density. The emission in the near 
infrared (NIR) is dominated by old populations - largely red giants 
- and will thus be an excellent tracer of the underlying mass 
distribution. However, some contamination arises from young 
stars - red supergiants (RSGs), and OB associations. Previous 
studies have shown these contributions to the flux to be limited 
to small spatial regions and even in these regions at most 20 
per cent of the flux comes from the 'contaminating' population 
JRix & Riekdl 19931) . In addition, the effects of interstellar dust are 
greatly reduced in the near infrared (studies of Ks band (2.2/im) 
data have demonstrated that the flux is attenuated by no more 
than 10 per cent, even in dust lanes dRix & Riekd 1 19930 ). The 
longer wavelengths detected by IRAC are likely to suffer even 
less attenuation. Despite the suitability of the near-infrared for this 
type of work, due to the inherent uncertainties outlined above it 
was decided to use several complimentary approaches to allow 
cross checks in our methodology. The primary data used for 



this research was o btained on the Spitzer Space Telescope b y 
the SINGS project JKennicutt et all l2003l : Iwillner et al.l |2004|) . 
IRAC bands 1, 2 and 4 (3.6, 4.5 and 8/im) are used; the 3.6 or 
4.5/im can be thought of as tracers of the underlying stellar mass 
distribution. Band 4, the 8 pva data, is used as a tracer of the shocks 
induced in the gas (further d iscussion in section |3). The IRAC 
data reduction is described in lRegan et al.l ( l2004h and the SINGS 
documentation. The IRAC data have a pixel scale (after drizzling) 
of 0.75 arcsec, and sensitivities 25.7 and 23.0 mag(AB) arcsec"^ 
(Icr) for 3.6 and 4.6/im respectively. The PSFs are approximately 
1.7 arcsec for 3.6 and 4.5/im, and 2 arcsec for the 8/im data. In 
addition, optical images in B, V, and / bands were used to create 
a mass map of the galaxy using a colour-correction technique 
described further in section 1231 The optical images were obtained 
on the Univ ersity of Arizona's B ok Telescope using the 90prime 
instrument (Williams et al. 2004j. Details on the data reduction 
can be found in Perez-Gonzalez et al. (2006). The pixel scale of 
the B, V, and / bands is 0.45 arcsec pix~^, with sensitivities of 
26.7, 26.3 and 27.9 mag(AB) arcsec"^ respectively. The PSFs 
are approximately 1.8 arcsec for all three band s. Finally, Ks band 
(2.2/im) images from the 2MASS catalogue ( IJarrett et all l2OO0h 
are used. The Ks band images have a sensitivity of 18.5 mag(AB) 
arcsec^^ and a pixel scale of 1.00 arcsec pix^^ (PSF ~3 arcsec). 
One arcsecond corresponds to ~17pc at 3.6Mpc, so at a radius of 
5kpc from the galactic centre, even with the largest scale images 
one pixel is equivalent to 0.1°. When compared with the error 
estimates for azimuthal angle quoted in section [3] it can be seen 
that the pixel scale is not a limiting factor in this analysis. 



2.2 IRAC 3.6 and 4.5^m data. 

To extract the stellar spiral structure from the 3.6^m data the ax- 
isymmetric components of the galaxy were first removed from the 
images. The non-axisymmetric components (residuals) can then 
be corrected for non-negligible small scale contamination by fore- 
ground stars and dust emission (features up to ~10 arcsec in di- 
ameter). The order of these steps is initially necessary because the 
correct adjustments require a visual inspection of the residuals (and 
so an axisymmetric model must be subtracted first). However an it- 
erative process would be possible, and in some cases preferable, in 
which the final fit to the axisymmetric components is carried out af- 
ter the small-scale contamination has been removed. In the case of 
M81 the difference is negligible because the star formation rate is 
relatively low, and the small-scale non-axisymmetric features make 
little difference to the fitting of the axisymmetric components. Al- 
though a stellar mass surface density map is not directly produced 
via the process described below, one can be easily created through 
the combination of the axisymmetric model with the fully cleaned 
residual images. 

T o fit the axisymmetric components, GALFIT dPeng et al.l 
120021) was used to produce a three-component fit to the galaxy us- 
ing a 2d implementation of a Sersic bulge (with index n=2.62 and 
Re=46.2 arcsec), exponential disc (Rs=155.4 arcsec) and constant 
background. The axis ratios (b/a) and position angles determined 
for the bulge were 0.71 and -31.9° and for the disc were 0.52 and 
-28.3° respectively. GALFIT has the ability to produce a model 
galaxy based on the best-fitting parameters, and this model was 
subtracted from the 3.6/im image to reveal the non-axisymmetric 
components, dominated by the spiral arms. The GALFIT model is 
discussed further in section[3] 

PAHs (polycyclic aromatic hydrocarbons) have emission 
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peaks, in additi on to continuum emission. A PAH emission fea - 
ture at 3.3^m ( JTokunaga etai]|l99ll : Ibulev & Williamslll98lh ) 
falls within the IRAC 3.6/im bandwidth, and so this is the most 
likely cause of much of the contamination in the residual image. 
PAH emission is also found in the S/im waveband, and this was 
used to reduce the contamination in the 3.6/im data by subtract- 
ing a scaled version of the 8/im data. The exact method is as fol- 
lows; the 8/im data were corrected for the stellar continuum emis- 
sion by the subtraction of a scaled version of the 3.6/im GALFIT 
model galaxy. The scaling constant used was 0.232, as given by 
iHelou et al.l J2OO4I) . It should be noted that, in using the GALFIT 
model (rather than the original 3.6/im image) to subtract the stellar 
continuum from the 8/im data, the continuum contribution from the 
spiral arms is not taken into account. However, the original image 
has been shown to contain PAH emission; by using the GALFIT 
model we remove the risk of affecting the PAH contribution to the 
8/im emission. The effect of this approximation (after all correc- 
tions), as can be seen in equation[T] is that the 3.6/im flux will be 
slightly lower on the spiral arms than if the full 3.6/im image been 
used to remove the continuum, but by less than a factor of ^^0.05. 
This small change is unlikely to be noticeable above the noise in 
the data. The systematic effect on the phase (if noticeable) will be 
to reduce the offset measured between the 8/im peaks and the den- 
sity maximum in the stellar spiral wave by slightly increasing the 
amplitude of the 8/im feature in phase with the stellar spiral. 

After correcting for the stellar continuum the 8/im im- 
age should have no remaining axisymmetric components. This 
continuum-corrected 8/im image can then be multiplied by a 
scaling constant and subtracted from the non-axisymmetric 3.6/im 
residuals to remove the PAH features. The scaling constant used 
was 0.08, chosen (by eye) to maximize the smoothness of the 
arms after the PAH subtraction, and this is judged to be accurate 
to ~ ±5 per cent. The process can be described by the following 
equation; 



^•^corr. — O.Drio 



n — ax2sym. 



OmiSdata ~ 0.232{3.6modsl)) d) 



The justification for the subtraction of PAH components from 
the 3.6/im image is relatively straightforward; a visual inspection 
of the data before and after the PAH features are removed shows a 
dramatic improvement in the smoothness of the arms (and any un- 
derlying spiral structure will produce a smooth mass distribution). 
The scaling constant used is in fact consistent with values suggested 
by models of PAH emis sion and observatio n s of interstellar emi s- 
sion the Galaxy (se e, e.g. lDraine &"Lil ( l2007h ; lFIagev et al.l ( l2006h ): 
iFlagev et al.l ( 120060 find a range in the ratio of 3.6/8/im intensities 
from 0.059-0.094 over 6 different regions observed, with an aver- 
age value of 0.072. 

It is likely that the flux from star forming regions also has 
small contributions from OB associations and RSG stars as dis- 
cussed in section 12.11 In addition, a number of foreground stars 
overlay the galaxy. The removal of remaining fine-structure was 
achieved by using the IRAF task XZAP, which removes any fea- 
tures which are na above the local background. The size of the fea- 
tures removed, and value of n, can be adjusted to give the most ap- 
propriate corrections. XZAP was slightly modified to use MEDIAN 
rather than FMEDIAN when calculating the smoothed image; the 
quantization of the data by FMEDIAN was found to have a notice- 
able effect in the fainter regions of the residual image. The best- 
fitting corrections turned out to have only a weak dependence on a, 
but the size of the smoothing box used was critical - too small and 
the small-scale features would be left virtually unchanged, whereas 



if the smoothing box was too large the features of interest could 
be affected. The best values were chosen after careful examination 
by eye of the final images in comparison to the original data. Final 
editing by hand was required for the largest foreground stars in the 
image, and this was carried out using IMEDIT in IRAF. Although 
the removal of stars in this way potentially leaves residual features 
from the PSF wings, these will be small in comparison to the PAH 
emission regions, and varying on scales that are much smaller than 
the spiral pattern. 

The IRAC 4.5/im image was also used, to provide a compari- 
son with the 3.6/im data. The 4.5/im data also appear to have PAH 
contamination, although with a smaller correction needed (a factor 
of 0.05 of the 8micron continuum-corrected image, as compared 
with 0.08 for 3. 6/im. This is consi stent with measured ratios of 
4.5/8/im flux in iFlagev et alj ( 120060 who find ratios in the range 
0.037-0.065, with an average of 0.048. Other steps in the image 
processing were exactly as for 3.6/tm and will not be described 
further. 

The process of creating a mass map can be seen in Fig. [TJ 
the top row shows the residual non-axisymmetric components af- 
ter removing the PAH components (for the IRAC data) and run- 
ning XZAP. These residual images are combined with the GALFIT 
models (not shown) to create the mass surface density maps on the 
bottom row. the original data is shown for comparison. 

2.3 Optical data. 

Optical B,V and / band images were used to create a stellar mass 
surface density map of the galaxy. The method relies on the rela- 
tionship bet ween mass-to-li g ht rat io (M/L) and optical colour as 
described bv lBell & de JongI ( l200ll) (BdJ) using B-V pixel-to-pixel 
colour to correct an / band luminosity map of the galaxy, using the 
relationship logio(M/Li) = a/ + bi(B-V). In their 2001 paper, BdJ 
give values for aj and bj as found from galaxy evolution models 
(the zeropoint, aj, assumes maximum disc M/L ratio). In this paper, 
as with PAH corrections for the IRAC data, the best-fitting value of 
bj is judged by eye so as to maximize the smoothness of the mass 
distribution. The value of b/ quoted by BdJ is higher than appears 
ideal in this case (see Table[T]l. The relative amplitude of the spiral 
structure (to axisymmetric components) turns out to be relatively 
insensitive to the value of b/ used; almost doubling the value of b/ 
from 0.275 to 0.5 only reduces the relative amplitude of the spiral 
arms by an average of ^^15 per cent, which is less than the dis- 
crepancies between relative amplitude estimates from the different 
wavelengths, as will be shown in section |42l In contrast, the value 
of aj only affects the zeropoint of the mass output, which is not im- 
portant when considering the relative amplitude of the spiral arms 
and so this was not investigated further. As a result, readers are ad- 
vised to exercise caution if creating normalised mass-models using 
the values of a/ and b/ found in this w ork. 

In a later paper. lBell et al.l ( 120030 offer a potential explanation 
for the discrepancy demonstrated in Table[T] with the acknowledg- 
ment that the original models do not use as large a metallicity scat- 
ter as observed in real galaxies. This metallicity scatter may serve 
to over-estimate the M/L slope, particularly in the NIR, although 
the corrections quoted do not reduce the M/L ratio in the / band 
sufficiently (the ratio decreases by a factor of ~ 1.3 rather than ~3 
as needed to reconcile the differences in Table [B. Ho wever, given 



the large scatter in the data (fig 20. iBell et al.N2003h ). the lower 
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M/L ratio found in this case would appear to lie within reasonable 
limits. 

The rationale for the choice of B, V and / bands is partly ne- 
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Figure 1. The non-axisymmetric residuals for (top, 1-r) IRAC 3.6 and A.S^im, Ks band, and / band colour-corrected mass maps respectively. The middle row 
shows the original data for the IRAC 3.6 and 4.6^m, Kg band, and the colour-corrected / band before the axisymmetric component is removed. The bottom 
row shows the mass surface density maps for each waveband in turn. 



cessity, since data were available in B, V, R, and / bands only. How- 
ever, BdJ recommend using an optical colour such as B-V or B-R, 
because these are effectively tracers of age (rather than metallicity 
as would be the case with a NIR colour). BdJ further note that a 
NIR band is best for mass estimates because the Ks or / bands have 
much smaller dynamic ranges than their optical counterparts. 

Once the mass map of the galaxy has been produced, the 
subsequent steps are as described above for the IRAC data; 
GALFIT is used to remove the axisymmetric component of the 



galaxy, and XZAP for star and fine- structure removal. 



2.4 2MASS data. 

The analysis of the 2MASS Ks band data was rather simpler than 
the previous two methods described above, possibly due to the 
relatively poor resolution and depth of the data (18.5 (AB)mag 
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Method 



a/ 



bj 



Bell & de Jong 2001 -0.627 1.075 

Bell et al 2003 -0.399 0.824 

Visual examination of these data -0.627 0.275 

Table 1. The comparison of the best-fitting values of the constants a/ and b/ 
as used in the method logio(M) = logio(Lj) -I- aj -I- bj(B-V). The value of 
a/ affects the normalisation of the stellar mass into solar units, and depends 
on whether the disc is sub-maxi mal or not. For this pape r the value of -0.627 
was chosen to be consistent with lBell & de Jong 120010 . 
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Figure 2. HI observation (known to trace shocked gas) on left, Spitzer S/^m 
on right, showing the degree of agreement between the emission at the two 
wavelengths. Images are to the same scale. The H I observation was obtained 
from NED, originally published bv iBraunI (1993). 



arcsec" for Kg compared to 25.7 (AB)mag arcsec" for 3.6/im). 
The two steps required were the fitting of an axisymmetric compo- 
nent by GALFIT and then XZAP to remove foreground stars and 
fine structure. 



2.5 IRAC S/iin data. 

The gas shock can be traced through the 8/im because this wave- 
length largely traces dust emission. A more conventional tracer of 
gas shocks is the HI 21cm line, but the 8/jm band has several advan- 
tages because any galaxy in the SINGS sample has high resolution 
8/im data available. Emission at 8/im is dominated by dust, which 
tends to be concentrated in regions of high gas density, which oc- 
curs at (or just behind) the shock front. To be visible in emission 
the dust needs to be heated, and the primary mechanism for con- 
centrated emission is the switch-on of young stars triggered by the 
shock front. The link between the two wavelengths is highly plau- 
sible (see Fig. ^, although not all 8/tm emission is from HI re- 
gions. The 8/im image is prepared for use by subtracting the stellar 
continuum emission, as described in section \T2\ but is otherwise 
unchanged. 



3 DISCUSSION OF MASS SURFACE DENSITY 
ESTIMATES. 

The two main quantities derived from the residual mass surface 
density maps are the phase and amplitude of the stellar spiral, and 
as such it is important to have a good grasp of the uncertainties in 
measuring these quantities. 

Over- or under-correction for PAH contamination via the 
scaled 8/tm image will lead to errors in the flux, and hence mass 
surface density. This is most likely to affect areas near the peak of 
the mass surface density because dust is concentrated in the shock. 
It is hard to quantify the possible error in the relative amplitude 
as a result of remnant PAH emission after all corrections, but it is 
possible to put a upper bound on the problem; a comparison of the 
relative amplitude for fully-corrected data with the measured rela- 
tive amplitude if no PAH or fine structure correction is used gives 
a difference of no more than ^^20 per cent, and at most radii no 
more than 15 per cent. The remaining error is likely to be a small 
fraction of this. The equivalent step for the optical images, colour- 
correcting the / band mass map will lead to an incorrect mass/light 
ratio and hence incorrect mass estimate. Dust lanes, star-forming 
regions and other features that are prominent in the B-V map will 
be particularly sensitive to the choice of M/L ratio, and as such this 
could introduce a systematic error in the peaks of the spiral arms 
and an incorrect value for the amplitude of the spiral. A scaling 
constant that is too large will over-subtract from the star forming 
regions, and increase the intensity in the dust lanes (the opposite 
is true if the scaling constant is under-estimated). In addition, the 
choice of smoothing box size in XZAP must be carefully consid- 
ered; if the smoothing box is made too large this has the potential 
to remove features that are part of the wave which would lead to 
systematic under-estimates of the amplitude of the spiral arms. In- 
vestigation into the effects of the choice of sky sigma and (pixel) 
width of radial samples in ELLIPSE demonstrates that the errors 
in relative amplitude are no larger than ~5 per cent for reason- 
able values. GALFIT was constrained to keep the same ellipticity 
and position angle (e and PA) values for all fits to all wavelengths, 
but scalelengths were initially allowed to vary. Fitted this way, the 
scalelengths for the two IRAC bands varied by < 1 per cent. In con- 
trast, the fit to the / band colour-corrected image gave a disc scale- 
length that was ~15 per cent smaller. When the fit is constrained 
to take the same scalelengths as the IRAC 3.6/tm fit, the relative 
amplitude of the spiral structure varies by no more than 5 per cent 
compared to the unconstrained fit over most radii, although the dif- 
ference is as large as ~15 per cent around R = 300arcsec. 

Examination of alternative fits to the PA showed that the sen- 
sitivity is ~ ±3°. Within this range the residual image is not no- 
ticeably affected, but the fit quickly appears worse as the PA is var- 
ied by larger amounts. As a further confirmation that the values 
used in this work are reliable, the PA of the exponential disc fitted 
by GALFIT agrees well with position angles dete rmined kinem ati- 
cally (e.g., within o ne sigma of the value found in lRotsI ( 119750 and 
within 2 degrees of lVisseil ( Il980bl) ). The fit in e is similarly good. 
As with the amplitude, constraining the / band scalelengths during 
fitting with GALFIT caused minimal changes to the phase; the fit 
to the gradient of the log spiral differed by no more than ~2 per 
cent over the region in which it follows a log spiral. 

As previously described, the PAH correction factor affects the 
mass distribution near the peaks of the spiral arms, which can create 
artificial trends in phase if the PAH correction is not correct. De- 
pending on the method used the phase can be quite insensitive to 
PAH contamination because these peaks are much sharper than the 
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underlying spiral. Using sine waves to fit the phase produces errors 
of no more than 7° between phase measured for fully corrected 
residuals and those with no PAH or small-scale structure coiTec- 
tions. The errors are not random, tending to steepen the gradient of 
the log spiral, thus reducing the pitch angle (by a maximum of ~5 
per cent in gradient). However, if instead a single point of highest 
intensity is used to measure the phase of the stellar spiral the error 
between phase estimates at each radius can be as much as ^^60° 
(discussed further in section |431 >. The / band colour correction has 
the potential to introduce systematic errors in phase because the 
offset of the dust lanes and stellar spiral varies with radius. How- 
ever, provided a suitable choice of M/L ratio is used this will not 
be a problem. The remaining small scale structures, mostly stars, 
are unlikely to bias the phase estimates. A small, well defined peak 
in intensity from a star is clearly distinct from the slow variation in 
underlying spiral structure. This is especially true where sine waves 
are being used to estimate the phase of the wave (see sectionQ, but 
even the methods of determining phase which are most susceptible 
to sharp peaks in intensity cannot introduce systematic errors into 
the results, since the positions of foreground stars are random. 

Finally, it is worth noting that if the 3.6 and 4.5/im images are 
compared, in addition to the large-scale similarities, the majority of 
small-scale (few pixel) structure is duplicated in both bands. The 
features observed are most likely individual AGB stars, meaning 
that much of the small-scale 'noise' is due to these stars. Since 
AGB stars bias the mass-light ratio, these stars will be the cause of 
much of the small-scale noise in the azimuthal profiles (see section 

n. 




Figure 3. Azimuthal profiles taken from the 3.6/im residuals. The range 
in azimuth is extended over Stt. Zero phase is along the -l-y axis in Fig.|4] 
increasing anticlockwise. The position angle of the disc is -28.3°. 



4 RESULTS. 

Given residual mass-maps with well defined errors it becomes pos- 
sible to analyse the data for density wave features. The IRAF task 
ELLIPSE was used to fit isophotal ellipses to a galaxy profile, and 
can then extract intensity profiles as a function of azimuth. The ax- 
isymmetric models of the galaxy produced by GALFIT were used 
to allow ELLIPSE to fit the position angle and ellipticity (PA and e) 
at different radial values (measured on the semi-major axis (SMA)). 
The e and PA values were then used to make elliptical intensity pro- 
files of the residual images. Ellipses were initially sampled every 

5 pixels along the SMA between 250-1000 pixels, with each pro- 
file averaging over the 5 pixels width to make use of all available 
data. These were then combined to create approximately logarith- 
mic steps, such that steps between successive profiles become 5, 
10, 15 or 20 pixels as radius increases. Finally, each profile was av- 
eraged azimuthally to give angular resolution ~ 1 ° . A sine wave of 
the form y — asin{2x + b) + c was then fit to each profile in order 
to extract information about the amplitude and phase of the stellar 
spiral structure. In addition to elliptical profiles of the galaxy, radial 
intensity profiles were obtained by sampling the residuals along 
narrow wedges (4° in angular size) from the centre outwards. 



4.1 Overview of trends In data. 

As can be seen in Figs.[3]and|4l the residuals are dominated by the 
central ring feature around 100 arcsec. The azimuthal profiles over 
the range 90<R<150 arcsec in Fig. [3] show approximately con- 
stant intensity except at the small breaks at either end of the ring. 
The profile at 150 arcsec also starts to detect the features that ap- 
pear to be short spiral arms (which can b e seen in Fig.[4t. The ring 
has been noted in previous work on M81 JGeorgiev & Getovll99ll : 



lElmegreen et al.lll995l) . but the two spiral arms that appear to con- 
nect to the ring have not been mentioned previously. The possibility 
that the inner spiral arms are purely an artefact due to subtracting 
the axisymmetric model (particularly errors due to the bulge fit) has 
been considered, but this does not seem to be the case: at 150 arc- 
sec the bulge contributes no more than a quarter of the light from 
the total axisymmetric components, and the ring and inner spiral 
can still be identified even if the bulge components are left in the 
residuals (and the disc alone is subtracted). The features are also 
immune to adjustments of a few degrees in the disc PA fits, beyond 
which the entire disc fit is obviously incorrect. In addition, there is 
sufficient agreement between the 3.6/im residuals and the structure 
in the 8/im (see, e.g. Fig.|4} to further support the identification of 
the inner spirals as real features. 

Between ~I50 and 300 arcsec the spiral arms appear to van- 
ish: although small perturbations with two-fold symmetry are visi- 
ble in this range the phase does not vary with radius; maxima are at 
~I30° - not aligned with the SMA (PA = 152°). The amplitude of 
the perturbations is small compared to the spiral arms beyond 300 
arcsec, with the relative amplitude increasing from 0.05-0.1. The 
cause of the variation in intensity is not clear, but does not appear 
to be a density wave. 

Beyond ~300 arcsec the spiral arms are prominent in the im- 
age, and the azimuthal profiles show clear sinusoidal oscillations 
with a phase dependence on radius. Thus the early indications sug- 
gest that a density wave is present in this radial range, and extends 
unbroken through almost 180°. 

In the 8/im image. Fig. [4] there is clear evidence of spiral 
structure. The grand-design spiral arms dominate outside ~300 
arcsec and can be traced over 180° of rotation. Inside this radius 
the structure becomes more filamentary and two-fold symmetry is 
less well defined. Relatively extended features are still visible, but 
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Figure 4. IRAC 3.6 and 8/im maps of non axisymmetric mass surface den- 
sity (residuals). Ellipses at 80, 150 300 and 675 arcseconds are marked. 
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Figure 7. IRAC 3.6/im, 4.5/im. and / band colour-corrected data showing 
the relative amplitude of the stellar spiral wave. The shaded area shows the 
error bars for the 3.6/.im data which are derived from the fits to the azimuthal 
profiles and are highly conservative. (Error bars for the 4.5/xm, and / band 
data are not shown, but are of a similar scale. 
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Figure 5. Azimuthal profiles taken from the 8^m image. The range in az- 
imuth is extended over 37r. Phase is defined as in Fig. [5] 



seldom seem to extend over more than ~90°. In addition to the 
large spiral patterns, there is a great deal of fine structure (spurs 
and feathering) visible in the 8/im image. These features can be ex- 
plained through t he shearing of small- s cale structures by differen - 
tial rotation (e.g IWada & Kodal ( l2004h : lOobbs & Boiineiil ( l2006h : 
iKim & Ostrikeil J2006l) ). 

There are some features that clearly show up in both the 3.6/im 
and 8/im wavelengths; the well defined spiral arms outside 300arc- 
sec are one such example, as are the much shorter spiral arms 
around ISOarcsec. There is a great deal of fine structure observed 
at 8/im that is not found in the stellar component. In the region be- 
tween 150 and 300 arcseconds it has been noted that, while the 8/im 
emission displays spiral-type features, there is little or no evidence 
for equivalent structure in the stellar distribution. The disjointed 
nature of the S/xm emission in this region may be a reflection of 
the fact that there is no strong stellar wave present. There is also 
no equivalent in the stellar mass distribution to the spurs and feath- 
ering observed in the gas (this is unsurprising given that the mass 
distribution is expected to be smooth over large scales). 



It is striking, although unsurprising, that the 8/im spiral fea- 
tures are much narrower than the corresponding features in the 
mass surface density maps, and are far less smooth when viewed 
along azimuthal profiles (see Figs. |3] and |5). As will be shown later 
(section l43t a sinusoidal wave can be used as a good approxima- 
tion to a density wave but this is clearly not the case for the gas 
response. 

Azimuthal profiles are more effective than radial profiles for 
studying the morphology of the stellar spiral structure; it is par- 
ticularly hard to determine the amplitude of the wave from radial 
profiles because the wave amplitude changes with radius. In con- 
trast, radial profiles are much more accurate for determining the 
position of the gas shocks from the 8/im data because the spiral 
is relatively tightly wound, and so the features are narrow when 
crossed radially, as can be seen by comparing Figs.|5]and|6] 

4.2 Amplitudes. 

The radial dependence of the amplitude of the stellar spiral struc- 
ture is shown in Fig. |7] The results show that both the IRAC and 
colour-corrected / band data display very similar trends of ampli- 
tude increasing with radius. The relative amplitude of the 2MASS 
data is probably suffering from noise effects; examining Fig. [8] 
shows that the errors in the radial range 300-600 arcsec are much 
greater for the 2MASS data than for the IRAC data. For this reason 
it was decided not to include the Ks band data in further analysis. In 
all data sets the relative amplitude is defined as (a/(model intensity 
+ c)) where a and c are the wave (half) amplitude and average in- 
tensity of the azimuthal profile respectively (as defined for the sine 
wave fit in section|4ll. The denominator is defined to correct for a 
non-perfect disc fit: whereas a perfect model would have averaged 
the disc intensity exactly to zero for all radii this has not happened 
(probably because the galaxy does not have a perfect exponential 
disc, or as a result of the spiral arms influencing the fit). 
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Figure 6. Radial profiles of the 8/jm image. The radius has been connected 
to face-on orientation in each case. 



4.3 Phase of the stellar spiral structure. 

The phase of the stellar spiral wave can be identified from the peak 
in the mass surface density on each spiral arm. The variation of this 
phase with radius reveals useful information about the pitch angle 
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Figure 8. IRAC 3.6/xm and Ks band data for the relative amplitude of 
the stellar spiral wave showing the sizes of the eiTors associated with each 
dataset. 








Figure 9. Estimates of the phase of the stellar spiral structure from peak 
finding (blue) sine wave fits (green) and estimates by eye (red) shown over- 
plotted on the 3.6/jm and S/iva images. It can be seen that the peak finding 
method detects contamination from spurs which are clearly visible at 8/xm. 



and behaviour of the spiral. The amplitude of the stellar wave is at 
most ~30 per cent of the mass of the disc and so the signal tends to 
get hidden in the exponential profile if the total mass surface den- 
sity is used. Hence, the easiest way to extract the phase information 
is to exclude the axisymmetric components by using the residual 
surface density maps. The sine waves fitted to the azimuthal pro- 
files, as described in section|4l provide one easy way of measuring 
the phase (if it is assumed that the spiral structure can be described 
by a perfect sine wave). The phase of the maximum is obtained au- 
tomatically from the fit, and a sample of fitted profiles are shown 
in Fig.[TT] As discussed in section |3] uncertainties still exist in the 
residual mass maps, and so a suite of complementary techniques 
been developed to identify the position of the potential minimum 
(or maximum in the mass) at each radius sampled. The method of 
picking the density maximum by fitting a sine wave was subject to 
a simple check by then examining each of the fits by eye, looking 
for approximately sinusoidal shape. For each profile the phases of 
the density maxima were identified (one point per arm). A final ap- 
proach was to use an automatic routine to pick the two highest peak 
positions (assuming the arms are ~180° apart). This method has 
the highest susceptibility of the three to bias from contamination 
by young bright objects (which will tend to create narrow, bright 
features). The automatic peak-finding method can be compared to 
the sine fits and fits by eye in Figs. l9l and 1 101 where it can be seen 
that, even after removing as much contamination as possible from 
the mass surface density maps, the peak finding method picks up 
traces of the spurs which are clearly seen in the 8/im image. 

In order to justify the fitting of the spiral pattern with just an 
m=2 sine wave, the power in other Fourier components was anal- 
ysed. Over the region of interest (~300-600 arcsec), the power in 
the next highest component, m=4, averaged only 6 per cent of the 
power in the m=2 components. Furthermore, examination of the 
fits demonstrated that much of this power was contributing to the 
features identified as remnant PAH emission: Fourier analysis of 
the non-PAH corrected azimuthal profiles showed that, on average, 
there was 40 per cent more power in the m=4 component than the 
PAH-corrected profiles (and ~30 per cent more in the m=6 com- 
ponent, the next highest increase). This can be seen qualitatively in 
Fig-EU where the effect of not removing the PAH component from 
the 3.6/im images is clearly shown. The figure shows profiles of in- 
tensity vs azimuth (plotted over a 37r range), with and without PAH 
corrections. It is clear that the deviations from a underlying sine 
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Figure 11. PAH corrected data (blue), and (red) XZAP but non-PAH cor- 
rected data at identical SMA values for the 3.6/^m wavelength data. The 
sine wave fits to the (PAH corrected) data are shown in black. The sharp 
peaks in excess of the sine wave are stronger in the non-PAH coirected data 
suggesting that the peaks seen in the PAH-corrected profiles are remnant 
PAH features. 



wave are much more pronounced in the non-PAH corrected data; it 
is likely that the remaining non-sine wave features in the azimuthal 
profiles are a result of not fully removing the PAH components. 

It can be seen from Fig.[TO]that although there is scatter in the 
phase, there is no systematic variation between the two arms, con- 
firming that the arms are symmetric and separated by 180°. The 
lower plot in Fig. [TO] illustrates the greater scatter in the phase 
determined from the peaks in intensity, particularly in the inner 
parts of the spiral (and is consistent with the presence of rem- 
nant PAH emission in the residual images). As already mentioned, 
this method has the greatest susceptibility to contamination, and 
although the data are shown for comparison they are not used to 
calculate the pitch angle or the offset from the peaks in the 8/im 
emission. It can also be seen that the phase only shows approxi- 
mately logarithmic spiral behaviour between 5.7^1n(R)^6.4, cor- 
responding to a range of 300^R^600 arcsec; the inner extent of 
the spiral wave must clearly lie at R ~300 arcsec. The 600 arcsec 
cutoff could mark the true extent of the spiral wave in the stel- 
lar population, but relatively low signal-to-noise by R ~700 arcsec 
might also play a part in the loss of log-spiral-like behaviour (see 
Fig- 13 ■ Using phase measurements from the sine wave and phase 
fits by eye the pitch angle can be determined, and is found to be 
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Figure 10. The top plot shows the two arms of the spiral wave with the tt phase difference removed, demonstrating that although there is scatter, there is no 
systematic variation between the phases of the stellar arms (they are 180° apart). The bottom plot distinguishes the sine wave fits and fits by eye from the peak 
positions found from maximum intensity, where it can be seen that there is much more scatter in the latter method. Zero phase lies along the SMA. 



23°. If the sine wave fits alone are considered (Fig|12l(. it can be 
seen that a steeper gradient (smaller pitch angle) can be fitted to the 
4.5/.im than to the 3.6/im or I band colour-corrected data. The devi- 
ations from the straight line (log spiral) behaviour are not insignif- 
icant; all three wavelengths show a tendency for phase to increase 
faster with radius at smaller radii than at the outer extent of the 
spiral. This causes deviations of 5 and 9° from the phase predicted 
by the best-fitting gradient for the 3.6^m at the inner and outer ex- 
tents of the spiral respectively (and similar for the colour-corrected 
I band), and 1 1° at the inner extent of the spiral for the 4.5/im data. 
The reason for such close agreement between two of the methods 
and not the third is not clear, but the difference in pitch angle ob- 
tained (~3°) is probably a good indicator of the uncertainties in the 
method. 



4.4 Phase of the S^m spiral structure. 

The position of the peaks in the 8/im emission as a function of ra- 
dius was identified initially by eye from radial cut data (as shown 
in Fig.[6]l. The data were plotted with phase as a function of ln(R), 
and a logarithmic spiral was assumed, allowing a best-fitting gra- 
dient to be determined. This spiral was then overlaid on the 8pm 
image, and (small) adjustments to the gradient were allowed until 
the match between the logarithmic spiral and the image was maxi- 
mized. The final fit has a pitch angle of 14°, shown in Fig. 1131 

The positions of the 8/im peaks can also be identified directly 
from the 8/im image, by using TVMARK in IRAF. The positions 
identified through this method are shown in Fig. [T3] The data are 
also shown in Fig.[T4l viewed as (j) vs ln(R). Interestingly the inner 
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Figure 12. The phase of the stellar spiral structure, as measured by fitting 
sine waves to the azimuthal profiles. The phase can be fitted with a logarith- 
mic spiral between 5.7^1n(R)^6.4, as shown. 



Lindblad resonance (discussed in section [J!2l ( seems to lie approx- 
imately at the radius of an abrupt phase shift in the data. 



4.5 Offsets between the stellar spiral and 8/^m emission. 

With well defined phases for the 8/im emission and stellar spiral 
as a function of radius it is possible to look for an offset between 
the two. When Figs. [9] and [T3] are compared it becomes apparent 
that there is a larger offset between the peaks of the best-fitting 
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Figure 13. Left; the S/im image with (blue) the positions of the peaks in 
emission as identified by eye from radial cuts, and the positions of the spiral 
fit to the peaks (red). Right; positions of the peaks in emission as identified 
by eye (red) and from the radial cuts (blue). 
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Figure 14. The positions of the peaks of the Sfiva spirals identified directly 
from the image (deprojected to face-on orientation) with the best-fitting log- 
arithmic spiral plotted for comparison. The vertical line at ln(R) ~5.5 marks 
the approximate position of the inner Lindblad resonance; although spiral 
structure in the gas can be fi'aced well inside this radius the density wave 
only exists at radii greater than the ILR. 



sine waves and the peaks in the 8/im emission than between the 
combined estimates for the density maximum (not exclusively esti- 
mated from the sine wave fits) and the 8/im. This is not surprising 
since, in section [43] it was shown that the sine wave fitting was 
the least susceptible to contamination from remnant PAH features 
which tend to be on or very near the peak of the 8^m emission. 

Fig. [T5] quantifies the angular offset for the different methods 
used to identify the phase of the stellar spiral and 8/xm emission. 
The offset is plotted using the offse t function, $, which for consis- 
tency with lGittins & Clarkd ( |2004|) (G&C) is defined as 2(08f,m - 
(p stellar -arm)- The calculatioH was performed using the data, not 
the logarithmic spiral fits; the only data excluded are those obtained 
from the (stellar spiral) phase peak finding method, which, as can 
be seen in Fig.|9l often detects spurs rather than the true spiral struc- 
ture. The trend is for the offset to increase with radius which is 



in agreement with predictions in G&C. The model closest to M81 
from G&C is plotted for comparison with this data; corotation is 
predicted to lie where the offset function reaches -n. 



5 DISCUSSION. 

5.1 Amplitudes. 

The radial dependenc e of arm am plitude has been examined in the 
past, (e.g. .Elmegreen et al.l ( Il989i) ). using the / band; these data are 
also plotted in Fig. [16] The results show that both the IRAC and 
colour-corrected / band data display a smoother trend with radius 
than the Elmegreen data, and the relative amplitude is generally 
smaller. This may be explained due to contamination effects from 
young stars in the Elmegreen data (there was no equivalent pro- 
cess to colour correction or PAH removal carried out). It is also 
worth noting that the relative amplitudes are defined differently; 
in this paper relative amplitude (in Fourier components) is F2/F0, 
theirs (F2-l-F6)/(Fo-l-F4-l-F8). However, components beyond m = 2 
are normally small so to first order these are unlikely to result in 
large differences. Other authors who have investigated the ampli- 
tude of grand design spiral waves include i Schweizen ([1979), who 
finds the relative amplitude of the spiral arms in the O band to vary 
between ~0. 1-0.5. For comparison, the relative amplitudes of the 
spiral structure in other grand design spirals have been determined; 
lElmegreen et al.l ( Il989i) include M51 and Ml 00 in their study and 
find that the relative amplitude in MlOO varies between 0.1-0.5. 
The arms in M5 1 appear to be somewhat stronger, varying between 
~0.2-0.8 (in both galaxies the overall trend is for the relative ampli- 
tude to increase with radius, as withMSl). Sch weizer finds the rela- 
tive amplitude in M5 1 to increase from 0.2-0.8. lRix & Riekd jl993h 
have also studied M5 1 with Ks band data and find a relative ampli- 
tude that varies between ~0. 1-0.5 , again obeying the gener a l trend 
to increase with radius. Finally, lElmegreen & Elmegreerj <199d) 
use optical data to constrain the spiral arms in NGC 1566, and find 
a density wave with an average relative amplitude of ~0.4. Values 
for relative amplitude used in simulation s of M8 1 have tended to be 
smaller than those found observationallv; ! Vissen 1 1980b.) used mod- 
els where the relative ampl itude ('forcing') varied between 0.05- 
0. 1 .i Gittins & Clarkd ( |2004|) use a model for M8 1 with relative am- 
plitude '^0.03, and investigate amplitudes no greater than 0.1 in the 
non-M81 -specific models. 

It should be noted that the large error bars associated with the 
relative amplitude (as shown in Fig. [7} are a conservative estimate 
of the error in measuring the amplitude of the wave. Given that 
the random scatter in the amplitudes for each wavelength is much 
smaller, the errors given should be considered to be an estimate of 
the systematic error associated with the measurement. 



5.2 Offsets. 



ICittins & Clarkd <2004) (G&C) used semi-analytic theory supple- 
mented by hydrodynamic simulations to compute the steady state 
response of isothermal gas to a rigidly rotating spiral mode. For 
given galaxy input parameters (i.e. rotation curve, relative ampli- 
tude of mode, azimuthal wave number) it is found that, inward of 
corotation, the shock in the gas moves steadily upstream with re- 
spect to the spiral arms in the stellar potential as radius is increased. 
As corotation is approached, the angular offset increases rapidly. In 
principle, therefore, offset data should readily indicate the corota- 
tion radius. In practice, however, the weakening of the shock as 
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Figure 15. The offset function between the 8/im emission (predicted to trace the gas shock) and density maximum in the stellar spiral as a function of radius 
(note, the offset function is twice the actual angulai' offset). Plots on the top line use the 8/im phase iden tified from radial cuts, whereas plots on the bottom use 
the 8/jm phase determined by eye from the Hfim image. In each plot the offset function prediction from lGittins & Clarke (2004) is shown as a solid black line; 
corotation occurs when the offset function reaches -7r/2, Note that the gas shock is predicted to lie upstream of the density wave as corotation is approached 
(as is observed). 



corotation is approached means that it is only feasible to obtain ob- 
servational constraints on offsets at radii well within corotation and 
G&C therefore present estimates of the accuracy with which the 
corotation radius can be determined as a function of the complete- 
ness of the data available. As a general guideline, G&C suggest that 
Rco can be determined to within 25 per cent if the angular offsets 
observed give an offset function range exceeding '^ 7r/4 (for a two 
armed spiral, offset function is twice the angular offset). 

If the observed spiral structure in M81 is assumed to be a long- 



lived density wave which maintains a well defined, constant pattern 
speed over a number of galactic rotations, it is possible to compare 
our results to those of G&C. When the data in Fig. [15] are exam- 
ined, corotation appears to lie at ~13kpc. This result has, for the 
moment, assumed that any discrepancies between the model pa- 
rameters in G&C compared to those determined in this paper can 
be ignored. This is largely justifiable due to small differences in 
the relevant parameters, however the difference between the rel- 
ative amplitude of the spiral perturbations in the model used by 
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Figure 16. Relative amplitude of the spiral structure determined for the 
IRAC an d / band data (as in Fi g. |7j compared to the relative amplitude 
found bv Elmegreen et alj 09891) . 



G&C and our data are rather larger. To compare the relative ampli- 
tude quoted in this work to the 'forcing' used to drive the simulated 
gas shocks in G&C it is n ecessary to inc l ude th e halo contribution 
to the axisymmetric mass. Ide Blok et al.l (2008) have used SINGS 
and THINGS (The HI Nearby Galaxies Survey) data to constrain 
the halo mass, and in the case of M81 find that the fraction of mass 
in the halo to that in the disc+bulge increases from ~0. 12-0.25 in 
the range 300^R^600 arcsec. Using this additional contribution 
to the axisymmetric components, the relative amplitude of the den- 
sity wave to total axisymmetric mass (equivalent to the forcing em- 
ployed by G&C) then varies between 0. 12-0.26. The standard value 
for the relative amplitude used by G&C was 0.05 at 8.5kpc. In con- 
trast, we find the relative amplitude to be closer to 0.225 at 8.5kpc 
Q As noted by G&C, for a given radius (and fixed corotation), in- 
creasing wave amplitude causes the offset to decrease rather than 
increase. Another factor that affects the offset function in the mod- 
els is the assumed sound speed in the gas: G&C found that at fixed 
R/RcH the shock moves further upstream in the case of warmer gas 
(see also iDobbs & Bonnelll ( 120070 who found that the shock was 
downstream even at relatively large R/Rcij for their very cold gas, 
although it should be noted that this result is also found in warm gas 
by G&C for the m =4 case - inapplicable to M81 - that is treated by 
iDobbs & Bonnelll) . The model of G&C plotted in Fig.|15|assumes 
a sound speed of 10 km s~^. We thus deduce that either the use of 
colder gas or the use of a larger spiral arm amplitude would cause 
the fraction R/Rcij at a given radius to be larger than predicted by 
the model shown in Fig. [TS] In addition, as shown in Fig. [7] the 
relative amplitude of the arms appears to increase with radius, and 
G&C note that a systematic variation in wave amplitude over the 
disc may lead to an inaccurate estimate for corotation. Taking all 
the uncertainties caused by relative amplitude into account, corota- 
tion is estimated at ~12kpc ±3kpc (estimating an error of 25 per 
cent according to the guidelines given by G&C). 

There have been many estimates made of the value of corota- 
tion in M81; a sample are listed in Table |2] The methods have var- 
ied, but brief descriptions of a number are offered below. IShu et alJ 



^ Note that there is some evidence from G&C that the shock moves further 
downstream at small radii in the case of larger density wave ampHtudes. 
This may explain the fact that the pronounced shift to positive offsets at the 
smallest radii in the M8 1 data is not reflected in the model. 



Author 


Rcfi (kpc) 


Elmegreen et al. (1989) 


9 


Gottesman & Weliachew (1975) 


11.7-12.6 


Lowe et al. (1994) 


10 


Roberts et al. (1975) 


11 


Rots (1975) 


11 


Sakhibov & Smirnov (1987) 


>12 


Shu etal. (1971) 


15 


Tamburro et al (2008) 


8.7 ±4.7 


Visser(1980) 


12 


Westpfahl (1998) 


9 


This work 


12 ±3 



Table 2. Values for Rcij found previously. The methods used are described 
further in the text. All values are scaled to the cun'ent M8 1 -distance estimate 
of 3.6Mpc. Errors are shown when stated by the authors, otherwise Rcfl is 
given to nearest kpc. 



( 1197 If ) uses the outermost H IT regions to define the radius of coro- 
tation. [Roberts et al.l ( 119751) use HII measurements from iMunchI 
(1 19591) in combination with spiral structure and the extent of the 
visible disc and, assume that Kcr must b e app r oxima tely co- 
incident with the radial extent of all three. iRotsI ( 119750 defines 
a pattern speed (f7p = 20kms~^kpc~^) which, when combined 
with the rotation curve, gives believable estimates for both the 
ILR and Rcfl- The estimates for the ILR and Rcfl assume that 
the density wave must be constrained between the ILR and Rcfl- 
HI morphol ogy is used to define th ese in ner and outer limits. 
In contrast, iGottesman & W eliachew ( 119751) use HI data to cal- 
culate the radial wavelength of the spiral pattern and, together 
with the e picyclic f r equen cy and mass surface density, use the 



witn tne e picyclic t r equen cy and mass surface uensity, use tne 
method of iLin et al. ( 119691) to calculate f2p and from that Rcfl- 

311 



ISakhibov & Smirnovl ( ll987D use HI velocity data. The axisymmet- 
ric rotation curve is modeled iteratively, and non-circular perturba 
tions are used to calculate Rcij (the uncertainty i n the result is due 
to the rotation curve being insufficiently defi ned) . 'We stpf ah J I' 1 99_ 
uses t he Tremaine & Weinbe rg method ((Tremaine & Weinbeg 
1 1984ft . lElmegreen et alj ( Il989ft use B and / band observations 
to identify resonance feat ures (specifically th e 4:1 resonance) in 
the spiral arm amplitu des. iLowe et alj ( 119941) also use data from 
lElmeg reen et al.' ('1989) to build a model of M81 based on modal 
theory. Tamburro et al. use the offset between the peaks of star for- 
mation and cold atomic hydrogen measured via 24^m and HI emis- 
sion respectively. The offsets are calculated for a number of radii by 
the cross-correlation of azimuthal profiles of the two wavelengths, 
and corotation is expected where this offset falls to zero. Finally, 
Ivissen ( Il980bl) chooses Q,p such that the ILR and Rcii have be- 
lievable values for the rotation curve used in his models (such that 
wave propagation is possible over a suitable radial range). 

Given the range of values quoted above, Rcfl = 12-13kpc is 
not inconsistent. Using this value, combined with the M81 rota- 
tion curve (see, e.g. Sofue et al. (1999)), the pattern speed is found 
to be -^nkms^^kpc"^. The inner Lindblad resonance (ILR) can 
also be estimated using the rotation curve and Q.p, and is found 
to lie at '^4.3kpc. The positions of corotation and the ILR can be 
seen on Fig.[T7] Corotation lies beyond the observable stellar wave 
and 8^m spiral structure, which is consistent with the theory that 
the shocks should weaken as corotation is approached. Th e ILR is 
expecte d to act as a barrier to density wave propagation dToomreJ 
( 1969l): lMarkl ( ll97ll) ). so if the ILR is indeed located at 4.3kpc the 
dramatic fall in amplitude at this radius can be explained by the 
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Name 



Rrad (kpc) Rciosest (kpc) Yr (kms 1) 



M82 
NGC 3077 



38.6 

48.4 



29.0 
18.6 



-250 
-50 




Figure 17. The ILR and corotation circles (12.8kpc and 4.3kpc) indicated 
on the 3.6/im residual mass surface density map image (left) and 8/im im- 
age (right). 



presence of the ILR, rather tha n modal amplitude dam ping as has 
been claimed in the past - e.g. lElmegreen et al.l ( Il989h . Although 
the ILR should mark the inner edge of (stellar) spiral structure, the 
inner mini-ring and stubby spiral arms at ~100 arcsec (described 
in section |4m i are well inside this radius. The most likely explana- 
tion is that these features are driven by the spiral structure in the 
gas (rather than vice versa), as the mini-ring and spiral arms clearly 
coincide with structures in the Smicron image - the ILR does not 
affect the gas, and so spiral structure is able to propagate through 
to the regions within. However, it is also possible that a weak oval 
distortion could provide the driving force. It should be noted that, 
as shown in figure [141 the spiral structure appears to become more 
filamentary for R<Rilr, which would seem to support these con- 
clusions. 

The agreement between the observed offset and theory is par- 
ticularly interesting because the models of G&C assume that the 
galaxy has had time to reach steady-state behaviour. If the spiral 
structure is highly transient the offset between density wave and gas 
shock would not have time to develop, hence, at least in the case of 
M81 it appears that the spiral arms must be reasonably long-lived, 
at least for the order of a few dynamical timescales. 



5.3 Fourier decomposition. 

The lack of power in the non m=2 components of the harmonic de- 
composition (described in section |43t is initiall y surprising, partic- 
ularly the relatively small contributions to m=4. lBertin et al.l ( ll989l) 
studied global modes in spiral galaxies and concluded that it is not 
easy to support single dominant modes of spiral structure in their 
models (whilst observing that such galaxies do exist), and suggest 
that two or three important modes are more stable. It is quite possi- 
ble that the external influences on M81 from M82 and NGC 3077 
(see section |5.4| l are the cause of the strength of the m=2 mode, 
although the excitation of modes by an external influence is not 
necessary fo r the existence of spiral structure in the modal picture. 
ISertin et al.l also noted that the stability of spiral structure over a 
number of rotation periods is reliant on the existence of one domi- 
nant mode; a larger number of modes will tend to cause the spiral to 
evolve gradually in a quasi-periodic manner. Thus, the observation 
of a single dominant mode in M81 is in fact consistent with the in- 



Table 3. Rrad is the projected distance to M81 from the companion galaxy. 
V jj is given relati ve to M81. Distances of closest approach are taken from 
the simulations of lYunl U999I) . The distance of the M8 1 group is taken to be 
3.6Mpc. 



ferred long-term stability of the spiral structure from the behaviour 
of the offsets between the stellar spiral and gas shock. 

5.4 Interactions and tidal driving mechanisms. 

Simulations bv lYunI ( 119991) based on HI observations demonstrate 
that M81, M82 and NGC 3077 underwent a three-way interac- 
tion approximately 2x10* years ago (nearest approaches of M82 
and NGC 3077 to M81 were 2.2 and 2.8x10* years ago respec- 
tively). The current (radial) velocities and projected radial distances 
of M82 and NGC 3077 with respect to M81 are given in Tabled 
along with the distances of nearest approach taken from Yun's sim- 
ulations. It is immediately apparent that NGC 3077 cannot be driv- 
ing the spiral structure in M81; the relative velocity is too small. 

In contrast, it seems entirely possible that M82 could provide 
the driving force behind the spiral structure observed in M8 1 . Given 
the limited information available all calculations must be approxi- 
mate, but if it is assumed that the pattern speed of the spiral must 
match the angular speed of the interaction between M82 and M81, 
the speed and radius of M82 are certainly of the correct order of 
magnitude; M82 is required to have a velocity relative to M81 of 
~500 kms~^ at perigalacticon to match the pattern speed of the 
spiral (the matching of pattern speed and angular interaction is not 
necessary, but the response of the stellar spiral will be stronger in 
this case). Further, the timescale of this interaction (a few xlO* 
years), matches the minimum predicted by G&C that is needed to 
observe an offset between the density wave and gas shock. 



6 CONCLUSIONS. 

We have used near infrared and optical images to produce mass 
surface density maps of M81 using IRAC 3.6 and 4.5/im data, op- 
tical B, V and / band images, and 2MASS Ks band data. We have 
shown that there is an underlying spiral wave in the old stellar pop- 
ulation of MSI. This spiral structure is not unbroken throughout 
the galaxy, but can be traced through nearly a full 180° rotation 
between 300<R<600 arcsec. The amplitude of the stellar wave is 
found to rise from '^O. 1-0.3 over this range. The amplitude esti- 
mates from all three methods agree to within the uncertainties, al- 
though systematic differences in the measured relative amplitude 
appear at larger radii. The pitch angle of the spiral can be deter- 
mined within this range and is found to be 23 ± 3°. 

By using the IRAC 8/im band as a indicator of the position 
of shocks in the gas we have measured the phase of the shocks in- 
duced by the stellar mass variations. The angular offset between 
the gas shocks and stellar spiral was used to determine the corota- 
tion ra dius, by using the method put forward by Gittins & Clarkd 
( 120041) . which assumes that the spiral pattern rotates as a rigid body 
(ie, a density wave). Corotation was determined to lie at 12 ± 3kpc. 
The position of the inner Lindblad resonance was extrapolated from 
the rotation curve and pattern speed, and found to lie at 4.3kpc. 
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These results are consistent with previous estimates, and with the 
observed morphology of the stellar spiral and gas structure. In ad- 
dition, the offset between the stellar spiral wave and gas shock sug- 
gests that steady-state behaviour must have been reached, meaning 
that a density wave with a constant pattern speed must have sur- 
vived for at least a few xlO* years. Further numerical simulations 
are required in order to investigate whether the observed offsets are 
consistent with the spiral structure originating from a galactic en- 
counter 

The dynamics and projected distances of MSl's nearest neigh- 
bours, M82 and NGC 3077, are examined, and by matching the 
pattern speed of the spiral arms with the angular speed of the com- 
panion galaxy (around M81), it is shown that the interaction be- 
tween M81 and M82 could potentially provide the driving force 
for the density wave. The timescale of this interaction (nearest ap- 
proach ~2.2xl0* years ago) is consistent with the (minimum) time 
needed to set up a steady state response to the density wave in the 
gas. 

A logical extension to this work is to extend the sample to 
more galaxies. There are many candidate spiral galaxies in the 
SINGS survey, and a larger sample will allow a larger range of 
questions to be addressed, such as investigating further the link be- 
tween stellar spirals, induce d shocks in the interstel lar gas, and star 
formation (as carried out bv lSeigar & JamesI ( 120021) ■ for example). 
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